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Abstract. We analyze damage nucleation and localization in the random fuse model with 
strong disorder using numerical simulations. In the initial stages of the fracture process, 
damage evolves in an uncorrected manner, resembling percolation. Subsequently, as the 
damage starts to accumulate, current enhancement at the tips of the microcracks leads 
eventually to catastrophic failure. We study this behavior quantifying the deviations from 
percolation and discussing alternative scaling laws for damage. The analysis of damage 
profiles confirms that localization occurs abruptly starting from an uniform damage landscape. 
Finally, we show that the cumulative damage distribution follows the normal distribution, 
suggesting that damage is uncorrected on large length scales. 



PACS numbers: 46.50.+a, 64.60.Ak 

1. Introduction 

Understanding the scaling properties of fracture in disordered media represents an intriguing 
theoretical problem with important implications for practical applications (1). Experiments 
have shown that in several materials under different loading conditions, the fracture surface 
is rough and can be described by self-affine scaling (0) with universal exponents (3). Scaling 
is also observed in acoustic emission experiments, where the distribution of pulses decays as 
a power law over several decades. Experimental observations have been reported for several 
materials such as wood (4), cellular glass (jsl), concrete (6) and paper (7), but universality 
does not seem to hold.The experimental observation of power law behavior suggests an 
interpretation in terms of critical phenomena and scaling theories, but a complete theoretical 
explanation has not been found. The statistical properties of fracture in disordered media are 
often studied with lattice models, describing the medium as a discrete set of elastic bonds 
with random failure thresholds (111). These numerical simulations are used in estimating the 
roughness of the fracture surface, which is found to be self-affine ^ [l^ jj,), and the 



power law distribution of avalanche precursors (12; J3,; J5;.16). While the results agree 
qualitatively with experiments a quantitative comparison is not always satisfactory. 

Apart from the comparison with experiments, an important theoretical issue is to 
understand the origin of the scaling behavior observed in the numerical simulations of lattice 
model. A very well studied model is the random fuse model (RE M), where a lattice of fuses 
with random threshold are subject to an increasing voltage (I2TI: [l|; Tl3l Ivh fisl: l20h . A 
resistor network represents a scalar analog of an elastic medium and is thus relatively simple 
to analyze, while retaining some important characteristic features of the problem. 

Simulations of the REM show that the type of behavior at macroscopic fracture is 
significantly influenced by the amount of disorder (flTh . When the disorder is narrowly 
distributed, materials breakdown without significant precursors. As the disorder increases. 



Percolation and localization in the random fuse model 



2 



substantial damage is accumulated prior to failure and the dynamics resembles percolation 
Indeed, in the limit of infinite disorder, the damage accumulation process can exactly 
be mapped onto a percolation problem ( 23). It is still debated, however, if percolation is 
applicable to the case of strong but non-infinite disorder. Recently, some evidence has been 
provided in this direction (24) suggesting that the critical exponent, u, of the correlation length 
in the RFM with strong disorder is same as that of uncorrelated percolation (i.e. u = 4/3), and 
that the fuse model is in the same universality class of percolation. Close to failure damage 
would then localize and the resulting crack roughness would ensue from a gradient percolation 
mechanism (24). 

In the present work, we propose a different interpretation of damage localization in 
the RFM: while in the initial stages damage accumulates as in a percolation process, the 
final crack nucleates abruptly due to current enhancement, yielding the observed localization 
profiles. As a consequence of this, the damage localization length follows closely the scaling 
of crack width. In addition, we test percolation scaling by simulating the RFM with triangular 
and diamond (square lattice with 45 degrees inclined bonds to the bus bars) lattice topologies 
fcll [ij) and different disorder distributions (uniform and power law). The numerical results 
allows us to exclude that u = 4/3, but could still be compatible with another value of u, 
although devations from scaling are observed for large sizes. 

Finally, we show that the cumulative damage probability distribution at failure can be 
collapsed for different lattice sizes and follows the normal distribution. This suggests that 
damage is the sum of uncorrelated variables and thus no long-range correlation are present 
as would be expected in critical phenomena. This fact, together with the abrupt localization 
occurring right at failure, suggests the validity of the first-order transition hypothesis discussed 
inRef. (13). 

The paper is organized as follows: in section II we define the model and in section III we 
discuss the role of disorder. The percolation hypothesis is tested in section IV, while section V 
is devoted to damage localization. The scaling of the number of broken bonds is discussed in 
section VI. Section VII presents the analysis of the failure probability distribution and section 
VIII summarizes briefly our conclusions. 

2. Model 

In the random thresholds fuse model, the lattice is initially fully intact with bonds having the 
same conductance, but the bond breaking thresholds, t, are randomly distributed based on a 
thresholds probability distribution, p(t). The burning of a fuse occurs irreversibly, whenever 
the electrical current in the fuse exceeds the breaking threshold current value, t, of the fuse. 
Periodic boundary conditions are imposed in the horizontal direction to simulate an infinite 
system and a constant voltage difference, V, is applied between the top and the bottom of 
lattice system bus bars. 

Numerically, a unit voltage difference, V = 1, is set between the bus bars and 
the Kirchhoff equations are solved to determine the current flowing in each of the fuses. 
Subsequently, for each fuse j, the ratio between the current ij and the breaking threshold tj 

is evaluated, and the bond jc having the largest value, max, 4^, is irreversibly removed (burnt). 
The current is redistributed instantaneously after a fuse is burnt implying that the current 
relaxation in the lattice system is much faster than the breaking of a fuse. Each time a fuse 
is burnt, it is necessary to re-calculate the current redistribution in the lattice to determine the 
subsequent breaking of a bond. The process of breaking of a bond, one at a time, is repeated 
until the lattice system falls apart. In this work, we consider a uniform probability distribution, 
which is constant between and 1, and a power law distribution p{t) oc t^^^^^, < t < 1 
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with/? = ^. 

Numerical simulation of fracture using large fuse networks is often hampered due to the 
high computational cost associated with solving a new large set of linear equations every time 
a new lattice bond is broken. The authors have developed rank- 1 sparse Cholesky factorization 
updating algorithm for simulating fracture using discrete lattice systems (20). In comparison 
with the Fourier accelerated iterative schemes used for modeling lattice breakdown (2S|), 
this algorithm significantly reduced the computational time required for solving large lattice 
systems. Using this numerical algorithm, we were able to investigate damage evolution in 
larger lattice systems (e.g., L = 1024), which to the authors knowledge, is so far the largest 
lattice system used in studying damage evolution using initially fully intact discrete lattice 
systems. However, in this paper, we consider results up to L = 512 due to insufficient 
number of available sample configurations for L = 1024, which will be considered in a future 
publication. 

Using the algorithm presented in we have performed numerical simulations on two- 
dimensional triangular and diamond (square lattice inclined at 45 degrees between the bus 
bars) lattice networks. For many lattice system sizes, the number of sample configurations, 
Nconfig, used are excessively large to reduce the statistical error in the numerical results (see 
Tables 1 and 2) . Each numerical simulation was performed on a single processor of Eagle (184 
nodes with four 375 MHz Power3-II processors) supercomputer at the Oak Ridge National 
Laboratory. The statistically independent Nconfig number of configurations were simulated 
simultaneously on number of processors available for computation. 



3. The role of disorder 



The disorder distribution has clearly an important effect on the fracture behavior and a 
classification has been proposed in terms of a scale-invariant spectrum (29; 30), which for 
the uniform thresholds distribution (0 < t < 1) in terms of intensive variables, at and ft (at), 
is given by dHHO) 

ftiat) = 2 -at for 0< at < 2 (1) 

where 

at = p (2) 

logL 

log LHpjt) 

Jt[at) = — : p (3) 

logL 

Hence, for the uniform distribution between and 1, the two control parameters 00, and 0oo 
that characterize the thresholds distribution p{t) close to zero and infinity, respectively, are 
given by 0o = 1, and 0oo = oo (29; 30), where 0o/oo are defined by 

. ito MW'A (4) 

t^O/oo y log(t) J 

For the power law thresholds distributions as the one we are using, the scale-invariant 
thresholds spectrum is given by 

ft{at) = 2- Pat where < < | (5) 



and the two control parameters are 0o = /5 and 0oo = oo (I29c l30r). According to Ref. (I29c l3( 
based on the values of these two control parameters, 0o and (poo, both the uniform thresholds 
distribution and power-law thresholds distributions (as long as /3 < 2 (30)) belong to the same 
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scaling regime characterized by diffusive damage and localization (see regime B of Figure 18 
in reference (30)). According to this analysis, if the exponent v were to be universal and is in 
the same universality class as that of uncorrected percolation as conjectured in (.24') . then we 
expect to find z/ = 4/3 for uniform and power law threshold distributions. 



4. Damage and percolation 

In the case of strong disorder, in the initial stages of damage evolution, bond breaking events 
occur in an uncorrected manner and thus resemble percolation. As the damage starts to 
accumulate, some degree of correlation can be expected because of the current enhancement 
present at the crack tips. A natural question to ask concerns the relevance of these correlations 
as failure is approached. If correlations are irrelevant one should observe percolation scaling 
up to failure, as in the case of infinite disorder. On the other hand, in the weak disorder 
case, the current enhancement at the crack tips is so strong that a spanning crack is nucleated 
soon after a few bonds (or even a single bond) are broken (,17.') . The interesting situation 
corresponds to the diffuse damage and localization regime, where a substantial amount of 
damage is accumulated prior to failure. In this regime one should test if damage follows 
percolation scaling up to failure. 

Hansen and Schmittbuhl have considered broad threshold distributions p{t) oc 
I'^+P^ < t < 1 with two different [3 = and 13 = \ values, to see if percolation scaling 
was observed for /3 > 0. Based on the similarities with percolation, they (l24l) suggested the 
following finite size scaling law for the fraction of broken bonds, given by 

Pf-Pc^L'^ (6) 

In the Eq. pj and Pc represent the fracture thresholds in a lattice system size of L and 
infinity, respectively. As the system size L ^ oo, the broken bonds at failure pf Pc- The 
correlation critical exponent u was found in Ref. (24) to be consistent with the percolation 
value z/ = 4/3. An additional test is provided by the damage standard deviation at failure A/ 
(|27|) which should scale as 

Af ~ L'^ (7) 



Here we test these scaling laws for a wider finite size range than Ref. (1241) . wherein 
simulations with sizes up to L = 60 are reported. The fraction of broken bonds for each 
of the lattice system sizes is obtained by dividing the number of broken bonds with the total 
number of bonds, Nei, present in the fully intact lattice system. For triangular lattice topology, 
Nei = (3L + 1)(L + 1), and for diamond lattice topology, Nei = 2L{L + 1). The lattice 
system sizes considered in this work are L = {8, 16, 24, 32, 64, 128, 256, 512}. However, 
since corrections to the scaling laws are strongest for small lattice systems, in the following, 
we use lattice sizes L > 16 for obtaining the scaling exponents. Table 1 presents mean 
and standard deviations in the broken bond density (fraction of broken bonds) at the peak 
load and at failure for various lattice system sizes in both the triangular and diamond lattice 
systems for uniform thresholds distribution. In order to compare diamond and triangular 
lattice topologies, we find it more convenient to use N^^i rather than L as a finite size parameter. 
This is because the two lattices have a different dependence of the real lattice size (i.e., N^i) 
on the linear size L. We plot in Fig. [H the mean fraction of broken bonds at failure as a 

function of N^i^^^ for diamond and triangular lattices, which in principle should obey Eq. |^ 
as well. While to accept the percolation hypothesis one should observe a linear regime, a net 
curvature is apparent in the data especially for sizes L > 100. A similar curvature is found in 
the 50% survival probability p^. 
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We repeated the same analysis using the power law threshold distribution with [5 = 1/20 
in triangular lattices of sizes L = 8, 16, 24, 32, 64, 128, 256, 512. Table 2 presents mean and 
standard deviations in the broken bond density (fraction of broken bonds) at the peak load and 
at failure for power law thresholds distribution. The result is reported in Fig. |2l and is once 
again in contrast with percolation scaling (i.e. u ^ 4/3) 

0.3 I • , • , • , 1 




Figure 1. The 50% survival probability ps (top) and the mean fraction of broken bonds pf 

(bottom) plotted as a function of N^i for the uniform threshold distribution. If percolation 
scaling is obeyed, the data should follow a straight hne. A net curvature is instead observed in 
all the data for large lattice sizes. 

Thus we can conclude that in the random fuse model with a uniform and power law 
distributions v is not equal to 4/3. On the basis of the results, however, we can not exclude 
the possibility that Eq. © is valid with a different value of u. This would correspond to some 
sort of correlated percolation, which, being a second order critical phenomenon, is expected 
to be universal with respect to the lattice structure. A way to test this idea is to plot the data 
as in Figs. [IE using a different value of u. In Figs.|3^ and we report similar plots for the 
triangular and diamond lattice topologies with uniform thresholds distribution, where it can 
be seen that a straight line is obtained only for large values of v. In addition, the v values for 
which a straight line is obtained in Figs.|3t and [Bj? are quite different from one another. For 
example, a direct fit of the data for L > 16 yields 

Pf - 0.0816 = 0.42 iVj°-i^ (8) 

for triangular lattices, and 

p/- 0.0751 = 0.57 iVj°-25. (9) 

for diamond lattices. Hence, Eqs. (HI) and ^ for triangular and diamond lattice topologies, 
estimate the scaling exponent u to be equal to 2.63 and 2.0, respectively. 
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Figure 2. The 50% survival probability and the mean fraction of broken bonds p f plotted 

as a function of N^i^^^ for the power law threshold distribution. Percolation scaling is not 
obeyed. 



An additional test is provided by plotting the standard deviation of bonds at failure, which 
in case of percolation should follow Eq.|7] For small lattice sizes the data follow a power law 
with an exp onent u ~ 1.56, but deviations occur at large sizes. In the framework proposed 
in Ref. (I24h one could interpret the data in Fig. |4l as an indication of an initial correlated 
percolation process with u = 1.56, which then crosses over to a localized state where scaling 
is obscured by the presence of a damage concentration profile. We will show in the next 
section, however, that there is apparently no localization at peak load (i.e. the maximum 
current before catastrophic failure). On the other hand, the damage standard deviation at 
peak load follows closely the behavior at failure as shown in the inset of Fig. HI Thus the 
damage profile does not appear to be responsible from the deviation from scaling. A different 
explanation would be that there is no scaling just because the lattice fails abruptly far from a 
(correlated) percolation critical point. 



5. Damage localization 

As discussed in the previous section, damage is not described by percolation critical scaling 
up to the largest sizes. Fracture is abrupt and damage localizes. Here we clarify when and 
how localization takes place. In particular, we will consider the damage accumulated up to 
the peak load (i.e. the highest current that the lattice can bear without breaking) and after the 
peak load up to failure. 

In Fig. 121 we display the damage profile p{y) at failure and at peak load in a single 
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Figure 3. The mean fraction of broken bonds pf plotted as a function of N^j^^^'^ for the 
triangular (top) and diamond (bottom) lattice with a uniform threshold distribution. A straight 
line is observed for different values of ly in the two case. 
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Figure 4. The standard deviation A/ of the fraction of broken bonds at failure plotted as 
a function of N^i for different lattices and disorder distributions. The curves are initially 
described by a power law with exponent —0.32 corresponding to v = 1.56, but flatten for 
L > 100. In the inset, we compare A/ with the same quantity computed at peak load in the 
case of a triangular lattice with uniform disorder. 



simulation, for the cases of uniform and power law disorders. For uniformly distributed 
thresholds, localization appears clearly at failure, while at peak load the damage spreads 
homogeneously though the lattice. For power law disorder, it is difficult to assess from 
these curves the extent of localization. To obtain a quantitative description of the localization 
process it is necessary to average over different realizations. 

Averaging the profiles is a delicate task since localization does not necessarily take place 
in the center, but can in principle occur anywhere along the length of the lattice. Thus one can 
not perform a simple average because this would yield a flat average profile irrespective of 
the individual profile shapes in a single realization. The authors of Ref. (26) proposed to first 
shift the profiles so that they are centered around the maximum and then average. This method 
emphasizes the noise too much, yielding a spurious cusp in the center. Another possibility is 
to shift instead by the center of mass of the damage or, to avoid any effects due to shifting, 
one can use the Fourier method. 

We first consider the damage accumulated up to peak load shifting the data by the center 
of mass method. The result displayed in Fig. |^ for the uniform disorder distribution clearly 
shows that there is no localization at peak load. In the case of power law disorder the profile 
as presented in Fig. IT] is not completely flat for small scales, but flattens more and more as 
the size is increased. We tend thus to attribute the apparent profile to size effects. These 
data imply that the localization profiles already observed in the literature (3; 25:l26l) reflect 
mostly the final breakdown event, after the final spanning crack is nucleated. To quantify 
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the corresponding damage localization, we can then average the profiles Ap{y) obtained 
considering only the damage accumulated between the peak load and failure. In addition, 
this procedure considerably reduces the background noise. 

The center of mass shifted averaged profiles for various system sizes are reported in 
Figs.EK and [8J5 for uniform and power law disorder cases and these profiles show that the 
profile shapes decay exponentially at large system sizes. The damage peak (Ap(0)) scales as 
^-0.3 j^Qj. uniform distribution and is roughly constant for the power law distribution. We 
have tried different ways to collapse these profiles to extract a localization length. The simple 
linear scaling {Ap{y, L))/(Ap(0)) = f{{y — L/2)/L), proposed in Ref. (24), does not yield 
a very good collapse (see Fig.|9K). A perfect collapse is instead obtained using the form 



where {Ap(0)) = and ^ ~ L", with a = 0.8 (see Fig.|9j5). The situation for power law 
disored is similar, but scaling is less precise. As shown in Fig.[8j5, the profiles are decaying 
again exponentially. We could not, however, perform a reliable collapse: a linear scaling 
seems appropriate for system sizes L < 256 but fails for larger sizes. 

To obtain additional confirmation of these results, we perform a Fourier analysis of 
the Ap profiles, thus avoiding any possible bias due to the shifting. We first compute the 
magnitude of the Fourier transform for each realization and then average over disorder. From 
Eq.[lO|we would expect the power spectrum of the profile to follow 



This result applies to an infinite system and finite size deviations and other problems of the 
discrete Fourier transform are expected to affect the data. Nonetheless, as shown in Fig.fTOl 
we can collapse reasonably well the curves using the same exponent as in real space. 

It is also possible to estimate the localization length directly, independently of the profile 
averaging. This is done by computing the width of the damage cloud asW = {{{yb—ybY)Y^'^, 
where y^ is the y coordinate of a broken bond and the average is taken over different 
realizations. We have first measured W at peak load and at failure, and obtained a result 
W ^ L, consistent with earlier results (29) (see Fig. fTTT). This result is expected, since 
for a uniformly distributed damage such as at the peak load, W ~ L/vT2 ~ 0.288L, 
and is in excellent agreement with the numerical data. This is consistent with the fact that 
there is no localization at peak-load and at failure the damage cloud is dominated by the 
uniform fluctuations already present at peak load. Any additional scaling due to localization is 
obscured. To uncover it we can restrict the average only to the bonds broken in the last failure 
event (post-peak damage), obtaining W ~ which is consistent with the data collapse of 
the profiles. It is interesting to notice that the width of the final crack scales as well with a 
similar exponent. This supports the idea that localization is produced by catastrophic failure 
which yields at the same time the crack and the damage profile. 

6. Scaling of the damage density 

It has been noted in the previous section that the final breakdown event is very different from 
the initial precursors. Thus, we consider the scaling of the number of broken bonds at the 
peak load, Up, that excludes the last catastrophic event. In Fig. [121 we plot Up as a function 
of the lattice size A'^p; for triangular and diamond lattices. The data displays a reasonable 
power law behavior Up ~ Nj^i, with b = 0.93 and b = 0.91 for triangular and diamond lattice, 
respectively as previously shown in Ref. dlSi) . Thus the difference between the two lattices 
is marginal and may be attributed to the results obtained from the smaller lattice sizes, where 
corrections to the fractal scaling may exist. By plotting Up/N^i vs N^i (see the inset of Fig. [T^ 



{Ap{y,L))/{Apm = f{\y - L/2\/0 



(10) 



{\p{k)\')/mm = f{kL'-'). 



(11) 
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Figure 5. Normalized damage profile p{y) — ^^j^^^^j^ at failure and at peak load in a typical 
triangular lattice simulation of size L = 512 for the cases of uniform and power law disorders. 
nb{y) denotes the number of broken bonds in the y*-^ section. 



we show that some systematic deviations appear. The data could be equally well fit by a linear 
law times a logarithmic correction Up ~ N^i/ \og{N(,i) as suggested in Ref. (26) (see Fig.fTBl). 
Both these fits imply that in the limit of large lattices the fraction of broken bonds prior to 
fracture vanishes (i.e. = in the thermodynamic limit). 

Figure [lUpresents the scaling for the number of broken bonds, {uf — Up), after crossing 
the peak load. Once again, the scaling exponent for (nj — Up) is similar for both the triangular 
and diamond lattice topologies, and is equal to 0.72 and 0.69 respectively. This exponent is 
consistent with the parameters estimated from the profile {uf — Up) ~< ]9(0) > (3L + 1)^ ~ 
L^-^, The behavior of power law disorder appears to be simpler. The number of broken bonds 
at peak load scales linearly with N^i, as shown in Fig. [15] A similar result holds asymptotically 
for {uf — Up) (see Fig.O, which is again consistent with the profile scahng. 

7. The failure probability distribution 

The cumulative probability distribution for the fraction of broken bonds at failure (also termed 
as cumulative failure probability distribution), defined as the probability 11/(^5, L) that a 
system of size L fails when the fraction of broken bonds equals = where nt, is the 
number of broken bonds, is plotted in Fig. [161 for varying triangular lattice system sizes. In 
Ref. (24) a data collapse of a different, but, related distribution (i.e., survival probability) was 
attempted using percolation scaling. As it is evident from the failure of such a scaling (i.e., 
Eq. (|6l)), this collapse is poor. 

On the other hand, we obtain a very good collapse by simply plotting the distribution 
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Figure 6. Average damage profiles at peak load obtained by first centering the data around the 
center of mass of the damage and then averging over different samples. The data correspond 
to unifrormly distributed disorder. 
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Figure 7. Average damage profiles at peak load obtained by first centering the data around the 
center of mass of the damage and then averging over different samples. The data correspond 
to power law distributed disorder. 




Figure 8. Average profiles for the damage accumulated between peak load and failure. The 
average has been performed after shifting by the center of mass. The profiles show exponential 
tails, (a) uniform disorder (b) power law disorder 




Figure 9. Data collapse of the profiles reported in Fig. [St- (a) Data collapse using a linear 
scaling for the localization length, (b) Data collapse using a power law scaling. 
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Figure 10. Collapse of the power spectra of damage profiles for uniform disorder. 
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Figure 11. The damage width at peak load and at failure are basically the same. The linear 
scaling is expected for a uniform distribution and it is not due to localization. On the other 
hand, localization can be observed for post-peak damage, the width scales as a power law 
similar to the one observed for the final crack. 
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Figure 12. Scaling of number of broken bonds at peak load for triangular and diamond lattices. 
The scaling exponents are very close to each other and the difference could be attributed to 
small size effects. There are, apparently, some systematic errors as shown in the inset. 
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Figure 13. The number of broken bonds at peak load can also be fit by a Unear function times 
a logarithmic correction by plotting n^jN^i as a function of N^x in a log-Unear plot. 
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Figure 14. The number of bonds broken in the last catastrophic event scales as a power law of 
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Figure 15. Scaling of number of broken bonds for a triangular lattice with power law threshold 
distribution. The scaling exponents for Up is very close to one and a similar result seems to 
hold asymptotically for — Up. 
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in terms of Vr = — = \ , where u„, and cr„, denote the mean and standard 

deviation of the number of broken bonds at failure, and pf and Ap^. denote the mean and 
standard deviation of fraction of broken bonds at failure. Fig. [l7] shows that Iif{p, L) may 
be expressed in a universal scaling form such that n/(p, L) = Ilf{pf) for both triangular and 
diamond lattice topologies of different system sizes L. A similar collapse can be performed for 
the power law disorder distribution (see Fig.fTSt. The excellent collapse of the data in Figs. \Tj\ 
and[T8] suggests that the cumulative failure probability distribution, Iif{pb, L) = Ilf(pf), may 
be universal in the sense that it is independent of lattice topology and disorder distribution.We 
have also checked that the distributions at peak load obey essentially the same laws, i.e., 
n(p) = lifipf) = Ilp{pp), where Iip{pp) is the probability that a system of size L is at the 
peak load when the fraction of broken bonds equals ph, and pp is the corresponding reduced 
variable at the peak load. Finally, the collapse of the data in Fig. fT9l indicates that a Gaussian 
distribution adequately describes 11/ . 

The fact that damage is Gaussian distributed suggests that there is no divergent 
correlation length at failure consistent with the conclusions of Ref. {!§) that reported a 
finite correlation length at failure. Long-range correlations in the damage would imply that 
the central limit theorem does not hold and hence the normal distribution would not be an 
adequate fit to the data. The absence of long-range corralation is again in agreement with the 
hypothesis that fracture is analogous to a first-order transition (jJJii). 




Figure 16. The cumulative probability distribution for the fraction of broken bonds at failure 
for triangular lattices of different system sizes. 
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Figure 17. The collapsed cumulative failure probability distributions for both triangular (t) 
and diamond (d) lattices of different system sizes with uniform disorder when plotted as a 
function of the reduced variable pf = {p — pf)/ Af. 



8. Discussion 

This paper presents numerical simulations on large two-dimensional triangular and diamond 
lattice fuse networks with uniform and power law disorder distributions. We focus our 
analysis on the statistical properties and localization features of the damage as a function 
of the lattice type, size and disorder form. The use of high statistical sampling and relatively 
large lattice sizes is essential to obtain reliable results. 

The picture emerging from our analysis is that, for strong disorder, damage accumulates 
first in the system in an uncorrelated (or short range correlated) manner. This process 
continues up to the peak load, where no apparent sign of localization is present. Further 
increase of the current leads to catastrophic failure through a large avalanche event, whose 
size scales with the lattice size as L^-'^^. The damage accumulated in this event is localized 
in a cloud surrounding the final crack. The damage profile has exponential tails and can be 
collapsed for different system sizes using an appropriate scaling law. 

The accumulated damage density at failure pf appears to be far from a percolation like 
critical point, since we could not find a reliable universal scaling law as the system size is 
varied. A different possibility that we have tested is that ^ as L ^ oo, but a fit using 
simple scaling forms show some deviation as well. Thus we can not conclude whether pf 
decays slowly to zero, or to a non vanishing asymptotic value. This value, however, would 
not necessarily coincide with a critical point. 

Finally, this study also presents the scaling of cumulative failure probability distribution, 
which is defined as the probability that a lattice system fails at a given fraction of broken 
bonds. Based on the numerical results presented, we show that the cumulative failure 
probability distribution is universal in the sense that it does not depend on the lattice 
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Figure 18. The collapsed cumulative failure probability distributions for lattices of different 
system sizes with power law disorder when plotted as a function of the reduced variable 
Pf — {p ~ pf)/Af. For comparison we also include two curves obtained in the case of 
uniform disorder. 



topology, i.e., the distributions are identical for both triangular and diamond lattice topologies. 
Furthermore, a normal distribution presents an adequate fit to the data. 
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Table 1. A summary of the main results of the simulations for uniform thresholds distribution, 
including the number of configurations used to average the results for each system size. Pp 
and p / denote the mean fraction of broken bonds in a lattice system of size L at the peak load 
and at failure, respectively. Similarly, Ap and A / denote the standard deviation of fraction of 
broken bonds at the peak load and at failure respectively. 
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Table 2. A summary of the main results of the simulations for the power law thresholds 
distribution, including the number of configurations used to average the results for each system 
size. Pp and p / denote the mean fraction of broken bonds in a lattice system of size L at the 
peak load and at failure, respectively. Sinailarly, Ap and A/ denote the standard deviation of 
fraction of broken bonds at the peak load and at failure respectively. 



